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With the dark energy phenomena explored over a decade, in this present work we discuss a specific 
case of the generalized Einstein-aether theories, in which the modified Friedmann Equation is similar 
to that in the Dvali-Gabadadze-Porrati (DGP) brane world model. We compute the joint statistic 
constraints on model parameters in this specific case by using the recent type la supernovae (SNe la) 
data, the Cosmic Microwave Background (CMB) shift parameter data, and the Baryonic Acoustic 
Oscillations (BAOs) data traced by the Sloan Digital Sky Survey (SDSS). Furthermore, we analyze 
other constrains from the observational Hubble parameter data (OHD). The comparison with the 
standard cosmological model (cosmological constant A Cold Dark Matter (ACDM) model) is clearly 
shown; also we comment on the interesting relation between the coupling constant M in this model 
and the special accelerate scale in the Modified Newtonian Dynamics (MOND) model initially given 
by Milgrom with the hope for interpreting the galaxy rotation curves without introducing mysterious 
dark matter. 



I. INTRODUCTION 

The independent discovery respectively in 1998 and 1999 that the current universe's expansion is actually speeding 
up rather than previously thought slowing-down due to the evolution dynamics dominated by cosmic matter compo- 
nent (mainly the speculated mysterious cold dark matter) [l| has been an amazing result. To account for that cosmic 
accelerating expansion, together with other astrophysics observations, such as CMB, large scale structure survey, like 
the SDSS, and the the universe age or Hubble constant measurements, a so coined dark-energy component (even 
more puzzling than the dark matter concept) with enough negative pressure has been hypothesized. According to the 
Wilkinson Microwave Anisotropy Probe (WMAP) 7-year data-set analysis it makes up about 72.8% of the uni- 
verse's contents. However, dark energy maybe the most mysterious component of the universe hitherto as envisioned, 
we know little about what it is and its nature. Over the past decade there have been many theoretical models for 
mimicking the dark energy behaviors, such as the simplest (just) cosmological constant and the popular quintessence 
models [3]. An alternatively instructive idea is that the general theory of relativity may fail to describe the universe's 
evolution on very large cosmic scales. Some attempts have been made at modifying the standard general relativity 
such as the f(R) extended gravity models string theory inspired cosmology models, brane cosmology and the 
holographic principle applicable to the universe's evolution modelings. In this present work we concentrate on the 
generalized Einstein-aether theories as proposed by T. G. Zlosnik, P. G. Ferreira and G. D. Starkman @, which 
is a generalization of the Einstein-aether theory developed by Ted Jacobson and David Mattingly 0, Q with a free 
general function T(K-) for model buildings. 

Arrangement for this paper is as follows. In the next section, Section UH we briefly review the framework of 
generalized Einstein-aether theory by providing its basic equations for our later use. In Section [ill! we take a specific 
form of the J-{JC) function allowed and discuss the corresponding modified Friedmann equations. In Section IIVI 
followed, we describe how to employ the observational data sets used for joint statistics analysis in details with the 
hope that this clear development can be useful to the related astrophysics and cosmology community. In Section [V] 
by figures and tables, we show our results compared with the currently standard cosmology ACDM model. The last 
section, Section [VII contains our conclusions and discussions. 
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II. GENERALIZED EINSTEIN-AETHER THEORIES 

In the early history for modern physics, the concept aether is considered to be a physical medium homogenously 
occupying every point in our universe. It determines a special rest reference frame, in which everything has absolute 
relative velocity respect to it. That suits for Newtonian dynamics very well. Later, this puzzling concept is rejected by 
Einstein's relativity with mainly optics experiments. By saying "aether" framework, a term in this present paper, we 
do not mean a mechanical medium naively, but rather a locally preferred state resting for each point of the spacetime 
in our physical evolutionary universe, determined by some hitherto unknown physics or its physical state is to be 
specified by some physical conditions with its environment. Some people even argue that the smoothly distributed 
CMB everywhere may be regarded as a modern version of aether. Einstein-aether theories were popularized by 
Gasperini in a series of papers A vector-tensor theory is suggested by Ted Jacobson and David Mattingly 0, Hf, 
where in addition to the metric tensor field of general relativity this theory also contains a time-like unit vector field 
which picks out a preferred frame at each point in the spacetime. Then it is generalized by T. G. Zlosnik, P. G. 
Ferreira and G. D. Starkman [!, H[. 

The action of this theory with the normal Einstein-Hilbcrt part action can be written in the form below 



S = Id X\J—g 



C-a + Cm 



(1) 



where La is the vector field Lagrange density while Cm denotes the Lagrange density for all other matter fields. The 
Lagrange density for the vector part consists of terms quadratic in the field and its derivatives [5[ : 

K, a \ a = c ig ^g 1(T + c 2 6ffi + c 3 5^, (2) 

where Cj are dimensionless constants and the coupling constant M has the dimension of mass. The FQC) is a free 
function that we do not know a priori, and the A is a Lagrange multiplier that enforces the unit constraint for the 



a/3 



7<r 



time-like vector field. In some papers an additional term c^A a A^ g ia is also included in the expression for K. 

We choose the inverse metric tensor g a @ and the contravariant vector field A 13 to be our dynamic degrees of freedom. 
Field equations from varying the action ([1]) with respect to g a P and A@ respectively are given by 

a/3 = J-a0 + Oirbrl a p (6) 
V a {F'j a p) = 2AA/3, (4) 

where T a p is the stress-energy tensor for the vector field and 

P = %> ja * = 2^ 7 VM 7 - (5) 

For the choice ([2]), T Q( g is given by [|[ 



1. 



1 



T Q(3 --V CT F'( J^A P) - J\ a A p) - J {aP) A°) -FY ioi0) + -g a pM 2 F+XA a A0, (6) 



2' 



where the sub (ah) means symmetric with respect to the indices involved and 

Yap = -c x [(VvAaWAp) - {V a A v )(y p A v )] . (7) 
In addition, the constraint that A is a time- like unit vector field gives A a A a = — 1. 

III. MODIFIED FRIEDMANN EQUATIONS 

Now we consider the case of a homogeneous and isotropic universe as preferred by the WMAP observations, which 
can be described by the Friedmann-Robertson- Walker metric 

ds 2 = -dt 2 + a 2 (t) ( l --^dr 2 + r 2 dQ 2 ) , (8) 

V 1 — kr z I 
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where k is the curvature parameter. In such a case, the vector must respect the spatial homogeneity and isotropy of 
the universe at large scales. Thus the only component non-vanishing is the time-like component. Using the constraint 
A a A a — —1, we can get 

A a = (1,0,0,0). (9) 
We take the matter component as a perfect fluid, so its energy-momentum tensor is of the form 

patter = ^jj^ + p (Jf aUp + g ^ ( 10 ) 

where U a is the fluid four- velocity. By using (jHJ and ((9]), JC can be simplified as 

where the coefficient a = c\ + 3c2 + C3 and the Hubble parameter H = a/a. It is easy to find by calculation that the 
stress-energy tensor for the vector field also takes the form of a perfect fluid, with an energy density given by (also 
see in f(J ) 

PA = 3aH 2 (F'-^) (12) 

and a pressure as 

PA =3aH 2 (-^' + I-)-apH -aT 1 -. (13) 
3 2/C a 

We can check that the vector field part's contributions obey the cosmological energy conservation relation pA + 
3H(pa + Pa) — 0. A simple case has been discussed by Sean M. Carroll and Eugene A. Lim in reference [lfj. Now 
we show that this conservation relation is applicable to an arbitrary form of .F(/C). 

By taking eqs.(l5))- pT|) into field equations © and ((U), the modified Friedmann equations can be derived (see also 
in I): 

1 aJ- . r t9 k 8nG . 
(l-aF' + - — )H 2 + - = — p 14 

2 /C a A 3 

^(-2H + aT'H) + ^r - 8irG(p + p). (15) 
at a A 

In order to see what has been modified, we list the standard Friedmann equations in the ACDM model below for 
comparison: 

tt2 k A 8vrG 

H 2 + — - T = -^~P 16 

a z 3 3 

-2^ + ^ = 8nG(p + p). (17) 
at a* 

We can see that a few terms involving T(1C) and its first order derivative are present, which may imply that an 
effective term involving cosmological "constant" or an effective cosmological "constant" (the effective vacuum energy 
for the universe) can be from the vector field's contributions, the mysterious "aether" . 

In their interesting works and [Hj |. a class of theories with .F(/C) — 7(— /C) n have been discussed. Noting that 
a ^ [12], and K is negative there. It has been shown that the scale M ~ H and for appropriate parameters the 
generalized Einstein-aether theories can lead to a late-time acceleration of the universe's expansion (for example, the 
n = case is just corresponding to the ACDM model). For that form of ^(/C), the first modified Friedmann equation 
(fT4l) becomes [|] 



H X 2(n-1)' 

M 



k 8nG 



H z + -~ = ~^P, (18) 



where e = — (2n — 1)7(— 2a)™/6. We can see if n — i, e = and the modified terms disappear. Thus, there will not 
be a modified term proportional to H . However, what about other forms of J-{K,)1 



4 



In the following part of this paper we consider a specific case, in which we take 

JY/C) = By^E+J— ln(-/C), (19) 
V a 

where (3 is a constant. Taking equation (fl"9]l into (fT4)) . Equation (|14l) then becomes 



fc 8ttG 



" - MH +- = —p. (20) 

CL O 

where we have used JC = 3aH 2 /M 2 as in given equation (fTTj) . For k = 0, equation (I2TJ1) is almost the same as that 
in DGP brane world model [l3l - il5l ]. For k ^ 0, it will be a little different. It is still unknown whether there is any 
relation between these two theories. 

It is easy to figure out that when p — > 0, H — > AI for the evolutionary universe whose geometry is almost flat at the 
late stage as indicated from the WMAP observations now. That is to say at late period of the universe's evolution 
when p oc a~ 3 — > 0, the universe will keep its accelerating expansion due to the existence of the aether field and an 
interesting scale appears M — > H then. 

Furthermore, we can calculate the effective state parameter for the vector field part and the deceleration parameter 
for our choice of the function .F(/C) : 

_ Pa H 

p A m 1 



ad H 

d = -Ta = -T72 - 1 = 3(«U + 1) - 1. (22) 

n £1 



From the (|22|) we know directly that to explain the speeding up of the universe's expansion as implied by lots of 
astrophysics observations, the effective state parameter for the vector field part's contributions today must be smaller 
than — | (instead of the — i as given directly from the standard Friedmann equations for the ACDM model). 



IV. CURRENT OBSERVATIONAL DATA 



A. Type la Supernovae 

The observations of Type la supernovae (SNe la) provide an excellent tool for probing the expansion history of the 
universe. Because all type la supernovae explode at about the same mass, their absolute magnitudes are considered to 
be all the same (M » —19.3 ±0.3). This makes them very useful as standard candles. The observation of supernovae 
measure essentially the apparent magnitude m. The theoretical distance modulus is defined as 



Mt/i = m th - M = 51og 10 D L (z) + /i , 



(23) 



where Dl(z) = Hq<1l{z) is the dimensionless luminosity and po = 42.38 — 51og 10 /i. Here h is the dimensionless 
Hubble parameter toady. Dl(z) is given by 



Dl(z) 



■ sinh 



H \/Q k 



dz' 
o W*) 



(24) 



where fifc is the fractional curvature density today. In this paper we use the Union2 data set consisting of 557 
supernovae. The corresponding Xsn function to be minimized is 



2 

Xsn 



557 

E 

i=l 



Hobs(Zi) — jJ,th{Zi 



(25) 



where 9 denotes the model parameters. The minimization with respect to po can be made trivially by expanding Xsn 
with respect to as [l(| 



XsN(0) = A-2p B + p 2 C, 



(26) 
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where 



557 



MO) = E 



i=l 
557 



s w = E 



i=l 
557 



Equation (|26[) has a minimum for /xq = B/C at 



X 2 SA r(0) = A(0) - 



B 2 (9) 
C ' 



(27) 
(28) 
(29) 

(30) 



Thus instead of minimizing Xsn we can minimize Xsn which is independent of /iq. 

B. Cosmic Microwave Background and Baryonic Acoustic Oscillations 

In addition to the type la supernovae data, we use the Cosmic Microwave Background (CMB) shift parameter and 
the Baryonic Acoustic Oscillations to compute the joint constraints. The shift parameter 7Z is defined in [ij} as 



where is the redshift of recombination and Da(z) is the proper angular diameter distance: 



D A (z) 



1 



: sinh 



dz' 



o H(z) 



(31) 



(32) 



The seven-year WMAP results [18| have updated the redshift of recombination = 1091.3 and the shift parameter 
K = 1.725 ± 0.018. The x 2 for CMB shift is 



Xcmb(0) 



[11(6) - 1.725]' 
0.018 2 



(33) 



Another constraint is from the Baryonic Acoustic Oscillations (BAOs) traced by the Sloan Digital Sky Survey 
(SDSS). In this paper we use only one node at z = 0.35. The distance parameter A is defined as[19j 



A = D V (Q.Z5) 



0.35 



where Dy is the effective distance 



D v (z) 



{l + zfD\{z) 



H(z) 



The value of A is given in [19]: A = 0.469 ± 0.017. Thus the X 2 for BAO is 

[A(0) - 0.469] 2 



Xbao(0) 



0.017 2 



(34) 



(35) 



(36) 



To compute the joint constraints, we add these \ 2 functions together: 



2 _ ~2 i 2 i 2 

X — Xsn + Xcmb + Xbao- 



(37) 
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C. Observational Hubble Parameter Data 



There are two major methods of independent observational H(z) measurement, which are called the "differential 
age method" and the "radial BAO size method". The details can be found in [20]. In that paper, Tong-Jie Zhang and 
Cong Ma summarize the up-to-date observational Hubble parameter data (OHD). See below in Table HI The data 
points at z = 0.24 and z = 0.43 are derived from the "radial BAO size method", while the others are derived from 
the "differential age method" as named. 

The x 2 for OHD is 

^ [H E( Zl ) - H obs (z,)? ( ^ 
»=o CT * 

where E(z) = H(z)/Hq is independent of Hq. Using the same trick mentioned before, the minimization with respect 
to Hq can be made trivially by expanding Xohd with respect to Hq as 

Xohd(0)=HIA-2H B + C, (39) 

where 

A - (40) 



1=1 

13 



B = ]P g (*) g jf( *) (41) 

i=l a * 

C = Y, H ° b j Zi) . (42) 

i—1 1 



Equation (f39|l has a minimum for Hq = B/A at 



xIhd = -^+C (43) 



Thus, instead of minimizing Xohd we can minimize Xohd which is independent of Hq. From Table HI we can see the 
errors of OHD (data sets) listed are relatively larger. So we do analysis only with the OHD (data sets) separately 
with the hope that we can obtain more accuracy OHD in the near future. 



TABLE I: The set of available observational H(z) data 



z 


H(z) ± lcr 
km s _1 Mpc _1 


References 


0.09 


69 ± 12 


[21, 22] 


0.17 


83 ±8 


[22] 


0.24 


79.69 ± 2.65 


[23] 


0.27 


77 ± 14 


[22] 


0.4 


95 ±17 


[22] 


0.43 


86.45 ± 3.68 


[23] 


0.48 


97 ±62 


[22] 


0.88 


90 ±40 


[22] 


0.9 


117 ±23 


[22] 


1.3 


168 ± 17 


[22] 


1.43 


177 ± 18 


[22] 


1.53 


140 ± 14 


[22] 


1.75 


202 ± 40 


[22] 
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V. RESULTS 



For the ACDM model, the general expression for the expansion relation can be directly written out as 



H(z) = H + fi fc (l + zf + fi m (l + z)3 + n R (l + z)\ (44) 
where Oa, fifc, ^m> are the fractional density of vacuum, curvature, matter and radiation today, respectively. 

" A =3iif' Qk= a>Hl> fim= W Qr= ^hJ ■ (45) 

In this paper, we fix = f2 7 (l + 0.227lN e ff) and take the present photon density parameter JX. = 2.469 x 10~ 5 h~ 2 
(for T cm b = 2.725-ftr) and the effective number of neutrino species at its standard value 3.04 [18|. We also use the 
prior h — 74.2 ± 3.6 given in [24j . So there are only two independent parameters fl\ and f2 m . The parameter Q k can 
be expressed by the others: 

n k = i - n A - n m - n R . (46) 

For the Einstein-aether theory with our choice of the free function J-(IC), we can solve H(z) from (|20l) : 



H(z) = H Q 

where we define 



v , +^(i + ^) 2 + fi m (i + ^) 3 + r! fi (i + ^ 4 



(47) 



M 

^ = tt-. (48) 



There are also only two parameters £Ia and Similarly, can be expressed as (|46[) . in which J7a needs to be 
replaced by Qa- 

Firstly, we compute the combined constraints from SNe la, CMB shift and BAO data sets. The results are shown 
in Table |TT]and Figure Q] The best- fit values for parameters of the ACDM model are consistent with the results given 
by WMAP 7 0. It gives a nearly flat universe geometry with a tiny fl k = 4 x 10 5 . For the Einstein-aether theory 
case as we choose, the best-fit f2 m is a little bit smaller than that in ACDM model, but it gives a larger fl k = 0.04. 
The ACDM model fits better to the data sets as its Xmin ^ s 26.182 smaller than that of Einstein-aether theory case 
as shown. 

Then, we do similar analysis with the OHD (data sets). The results are shown in Table HT1 and Figure [Has well. 
This time, we can see that the Einstein-aether theory case fits a little better to the data sets, and the difference of 
the Xmin i s ver y small now. 

Noting the similarity of our modified Friedmann equation to that in DGP brane world model and data analysis 
being done to DGP model [111, [26[, these result are rather natural, because best-fit fl k is very small. 



TABLE II: Best-fit parameters for ACDM model and Einstein-aether theory 





Model 


Best-fit parameters 


/Cmin 


SN Ia-CMB shift-BAO 


ACDM 


n m = 0.272 


542.693 






n A = 0.728 






Einstein-aether theory 


n m = 0.220 


568.875 






Qa = 0.739 




OHD 


ACDM 
Einstein-aether theory 


n m = 0.322 
fi A = 0.807 
VL m = 0.301 
Q.A = 1.040 


8.046 
7.989 



Using the results of combined analysis we also plot the effective state parameter of the vector field part's contribution 
wa{z) (Figure [3]) and the corresponding deceleration parameter q(z) (Figured]), for demonstration. 




FIG. 1: Joint constraints from SNe la, CMB shift and BAO. The la, 2a, 3a confidence interval contours of fi m and Qa (or Ha 
for the Einstein-aether theory case) in ACDM model (left) and the Einstein-aether theory case(right). 




FIG. 2: Constraints from OHD. The la, 2a, 3a confidence interval contours of Q m and Ha (or Qa for Einstein-aether theory 
case) in ACDM model (left) and Einstein-aether theory case (right). 




FIG. 3: The effective state parameter of the vector field part's contribution. 
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FIG. 4: The deceleration parameter. The thick solid line is the result of the Einstein-aether theory, while the thin dashed line 
represents ACDM model. 



VI. CONCLUSIONS AND DISCUSSIONS 

In this paper we only consider in details a specific case of generalized Einstein-aether theories and compute the joint 
constraints from observations such as SNe la, CMB shift, BAO data sets, and OHD respectively. Even though we only 
investigate a specific case of generalized Einstein-aether model, we already see that it has shown lots of interesting 
features, by well fitting to the combined data sets of the SNe la, CMB shift and BAO, as well as OHD respectively, 
and comparing with the ACDM model. 

The observational Hubble parameter data we have possessed now are relatively few and not so accurate. However, 
with the improving quality of observational H(z) data and more data points being measured (more sample compiled 
hopefully), it will be certainly a directly useful tool to test dark energy models and modified theories of general 
relativity, as well as corresponding cosmology models. 

For the case we analyze in this paper, the modified Friedmann equation is similar to that in DGP brane world 
model. It may be caused by the special T{K.) we choose. It looks first rather strange if there is any possible relation 
between these two theories, because the Einstein-aether theory is initially proposed for possible Lorentz violation and 
preferred frame effects, while the DGP brane world model considers a 3-brane embedded in a 5D bulk space-time. 
However, it is not completely impossible now. 

Moreover, it is clearly shown in this special model that M = HaHq ~ Hq, which is consistent with the requirements 
of MOND limit 0, HI, [Tl|. However, further work needs to be elaborated on the stability analysis of this specific case 
we have considered. 

For the last point of this present work we would like to make (but not the least importance) , we should emphasize 
especially that theoretically we can not give for granted that the phenomenological MOND theory can reproduce all 
the systematics of Rotational Curves (RCs) observations, although the MOND model fits BETTER than the ACDM 
based mass models. It is still far away to reproduce the wide and far telling systematics of the spiral galaxies' RCs 
(and of the mass distribution in the corresponding galaxies) [27J , so there will be lots of detail work to be done with 
any modified gravity proposals at least on the galaxy scale. 
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